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Abstract 

The newly operational X-ray satellite Suzaku observed the supernova remnant (SNR) RCW 86 in 
February 2006 to study the nature of the 6.4 keV emission line first detected with the Advanced Satellite 
for Cosmology and Astronomy (ASCA). The new data confirms the existence of the line, localizing it for 
the first time inside a low temperature emission region and not at the locus of the continuum hard X-ray 
emission. We also report the first detection of a 7.1 keV line that we interpret as the K/3 emission from 
neutral or low-ionized iron. The Fe-K line features are consistent with a non-equilibrium plasma of Fe-rich 
ejecta with n e t < 10 9 cm -3 s and kT e > 1 keV. We found a sign that Fe K a line is intrinsically broadened 
47 (35-57) eV (99% error region). Cr-K line is also marginally detected, which is supporting the ejecta 
origin for the Fe-K line. By showing that the hard continuum above 3 keV has different spatial distribution 
from the Fe-K line, we confirmed it to be synchrotron X-ray emission. 
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1. Introduction 

Since synchrotron X-ray emission was for the first time 
localized to a shell of a supernova remnant, SN1006 

(Koyama et al. 1995), a number of young SNRs have 
been proven to show synchrotron X-rays from their shells 

[see for example, Cas A (Hughes et al. 2000); Tycho 

(Hwang et al. 2002); or G347.3-0.5 (Slane et al. 2001)]. 
Synchrotron X-ray emission is from electrons with ener- 
gies up to 10 TeV and therefore, is the strongest evidence 
to date of particle acceleration in the diffuse shocks of 
young SNRs. 

RCW 86 is one of the SNRs from which hard X-ray 
emission, which can be synchrotron emission, has been 

discovered. It has, however, a distinct characteristic to 

have hard X-ray emission from a south-west part of the 
remnant associated with a neutral Fe-K line emission. 
Since the line is discovered in 1997 in the ASCA observa- 
tions of the remnant (Vink et al. 1997), its origin remains 

a mystery and the origin of the hard X-ray emission has 
also been unknown: synchrotron or another mechanism. 


In the initial paper, Vink et al. (1997) tried to explain pe- 
culiar combination of abundances, the strong 6.4 keV-line 
and weak line emission of light elements, by an electron 
distribution with a supra-thermal tail. BeppoSAX ob- 
servations confirmed the existence of 6.4 keV line from 
the south-west region of RCW 86 (Bocchino et al. 2000). 
They found spatial variation of hardness in the SW region, 
and explained the peculiar abundances by mixture of two- 
temperature plasmas: low abundance low-temperature 
and high abundance high-temperature plasmas. They at- 
tributed the latter plasma to Fe-rich ejecta. The detec- 
tion of high-energy (> 10 keV) emission by RXTE (Petre 
et al. 1999) suggested the existence of non-thermal X- 
ray emission. Bamba et al. (2000) and Borkowski et al. 
(2001a) found that the ASCA spectrum from the SW re- 
gion can be explained by a combination of three com- 
ponents: low and high temperature plasmas, and non- 
thermal (synchrotron) emission. Using the superior high 
spatial resolution of Chandra, Rho et al. (2002) spatially 
resolved the low-temperature plasma and the hard contin- 
uum emission. By showing spatial correlation between the 



hard continuum and radio emission, they confirmed that 
the hard continuum originates from synchrotron emission. 
Since they found the Fe Ka line only from the region of 
the hard continuum, they suggested that the synchrotron 
X-ray is emitted from Fe-rich ejecta. Hard X-ray emission 
from a north-east region of RCW 86 shows no evidence for 
a 6.4 keV line, and Vink et al. (2006) recently reported 
discussions on the ground that the origin of that emission 
is synchrotron. 

Since spatial distribution of the 6.4 keV-line and precise 
spectroscopy of the line features will give us important 
information on the emission mechanism of the line, we 
performed a Suzaku observation on the south-west part of 
RCW 86. In this paper, we mainly report on the results 
in the hard band (> 3 keV); hard continuum, and Fe- 
line features. We assume a kinematic distance of 3 kpc 
(Rosado et al. 1996). 

2. SUZAKU Observation of RCW 86 

Suzaku is a joined Japanese-US mission that was 
launched in July 2005 (Mitsuda et al. 2006). Two in- 
struments, the XIS and HXD, are operational. The XIS 
consists of 4 independent CCDs cameras, three front- 
illuminated (energy range 0.4-12 keV) and one back- 
illuminated (energy range 0.2-12 keV). Each is located in 
the focal plane of a dedicated nested thin- foil X-ray tele- 
scope (XRTs). The total effective area is quite large (~ 
1000 cm -2 at 6 keV). In addition, the low-Earth orbit of 
Suzaku contributes to a low and reproducible background. 
The HXD is a non-imaging collimated detector which ex- 
tends the bandpass of the observatory up to 600 keV. More 
information on the mission and the instruments on board 
(XRT, XIS, and HXD) can be found in Serlemitsos et al. 
(2006), Koyamaet al. (2006), and Takahashi et al. (2006). 

The south-west region of RCW 86 was observed by 
Suzaku in February 2006. At the writing of this paper, the 
background determination for the HXD is still not com- 
plete and we chose not to include any HXD data in our 
analysis. The XIS contamination (Koyama et al. 2006) 
plays no role at the energies of the Fe-K line, our main 
concern for this paper. 

The data (revision 0.7; Mitsuda et al. 2006) were an- 
alyzed after the following filtering. Data taken with low 
data rate, at low elevation angle from the Earth rim (<5°), 
or during passage through the South Atlantic Anomaly 
were removed. We also removed data from hot and flick- 
ering pixels using the cleansis software. In order to sup- 
press the hard-band background of XIS 1, which is higher 
than that of XIS 0, 2, or 3, data taken at low cut-off rigid- 
ity regions (< 6 GV) were also removed from the XIS 1 
data additionally. The resultant exposure times are 116 ks 
for each XIS 0, 2, 3 and 96 ks for XIS 1. At this stage 
of the mission, the attitude is not completely cbrrect and 
we used a point source detected in an XMM-Newton ob- 
servation of RCW 86, to correct the absolute coordinate 
system of our observation. The correction is rather small 
(ARA=— 0.°0044, ADec=— 0.°0049) for all 4 of the XIS. 


3, XIS Images 

Figures la and b show XIS intensity contours in the 
low (0.5-1.0 keV) and high (3. 0-6.0 keV) energy bands. 
The high-energy band was chosen to show the contribu- 
tion from the continuum only. We have shown on Fig. la, 
the 6.4 keV-line intensity image. This image was gener- 
ated by subtracting from the (6.34-6.46 keV) image, the 
continuum contribution extracted from the (5.0-6. 2 keV) 
image. Relative normalization between these images is 
done using a power- law model. Since smoothing may show 
artificial structures, we adaptively binned the image using 
the weighted Voronoi tessellation algorithm implemented 
by Diehl & Statler (2006), which is based on the algo- 
rithm of Cappellari & Copin (2003). The size of each bin 
was determined so that significance of the 6.4 keV line is 
at least 5 a. The resulting image is shown in Fig. lb. In 
order to quantify spatial correlation between the 6.4 keV 
line and the continuum emission, we furthermore obtained 
spatial variation of count-rate ratios between the 6.4 keV 
line and the (5. 0-6.2 keV) band using the same binning. 
The result is shown in Fig. 2 with the 3.0-6. 0 keV intensity 
contours. If we multiply the ratio by 3, it corresponds to 
the equivalent width (keV) of the 6.4 keV line at that po- 
sition. We have generated in Fig. 3 a true color image of 
the remnant showing the spatial variations between low 
and high-energy continuum (shown in blue and red re- 
spectively), and the region of emission of the 6.4 keV-line 
(shown in green). 

The spatial distribution of the 6.4 keV-line is revealed 
for the first time thanks to combination of the large ef- 
fective area and the low-background of Suzaku at the en- 
ergy band around 6.4 keV. The red region in Fig. 3 cor- 
responds to the location of the low-temperature plasma, 
which was attributed to the forward shock by Rho et al. 
(2002), whereas the blue region corresponds to the fila- 
mentary structures of the hard X-ray emission revealed by 
the Chandra observation (Rho et al. 2002). The 6.4 keV- 
line is emitted from the inner part of the red region, or 
the forward shock. As apparent in Fig lb, the high-energy 
continuum emission from the remnant does not correlate 
with the position of the maximum intensity of the 6.4 keV- 
line. Fig. 2 shows that the equivalent width of the 6.4 keV 
line and the continuum emission rather anti-correlate with 
each other, the origins of the line and the continuum are 
strongly suggested to be different. We will examine be- 
low how the spectral analysis of the line fit with this new 
finding. 

4. Spectral Analysis of the 6.4 keV Line 

In order to study the origin of the 6.4 keV emission 
line, we extracted spectra from a region shown in Fig. 1, 
an elliptical region of 12.'5 x 9.'2 diameters with its center 
on (RA, DEC)j 2 ooo == (14 h 41 m 08?7,-62°40'10") inclined 
for 20°. This region was chosen so that maximum signal- 
to-noise ratio can be obtained for the 6.4 keV line. Even 
though one may expect three distinct spectra from three 
color regions in Fig. 2, the spatial resolution of XIS did 



Table 1. Suzaku Observations Appearing in this Paper. 


Target 

Observation ID 

Date 

Exposure 3 - (ks) 

RCW 86 Southwest 
North Ecliptic Polar 
High Latitude A 
SgrC 

500004010 

500026010 

500027010 

500018010 

2006/02/12-14 

2006/02/10-12 

2006/02/14-15 

2006/02/20-23 

116.5/96.4 

96.0/87.5 

79.5/71.1 

114.5/98.6 


Exposure times for XIS0,2,3/XIS1 after the data screening. 



Fig. X. XIS 0.5-1 .0 keV (a) and 3. 0-6.0 keV (b) intensity contours overlaid by the 6.4 keV-line images. The 6.4 keV 
line image in (a) is smoothed with a Gaussian kernel of cr=25", whereas that in (b) is adaptively binned using the 
weighted Voronoi tessellation algorithm (Diehl & Statler 2006). Images and contours are in linear scale. The spectral re- 
gion is shown with a broken-line ellipse. The XIS field of view is designated with a solid-line square in each image. 


not allow us such distinct separation. Moreover Chandra 
and XMM-Newton are better in spatially resolved spec- 
troscopy for the soft component and the hard continuum 
emission. Therefore we here concentrate on the spec- 
troscopy of the Fe-K line features. 

As the entire field of view is covered by the remnant, 
we used ” blank-sky” observations, extracting the back- 
ground from the same regions on the detectors as the ones 
used in the analysis. Both observations (’’North Ecliptic 
Polar” and ’’High Latitude A”, see Table 1) were carried 
out either immediately before or after our observation of 
RCW 86, insuring an identical response of the detector. 
Spectra of these two observations were weighted averaged 
using the exposure time. Above 3 keV, differences be- 
tween the responses from all the FI chips are negligible 
and we summed the spectra of XIS 0, 2 and 3 for the 
following analysis. Although the BI is kept and analyzed 
separately, the differences between the summed FI and 
the Bl-derived spectral parameters are negligible. 

Background-subtracted spectra in the energy band 3.0- 
10 keV are shown in Fig. 4. Since we see power-law like 
spectrum with two conspicuous lines at 6.4 and 7.1 keV, 
we fitted the spectrum with a power-law function plus 
two Gaussians corrected for the interstellar absorption. 


Since absorption column is not well determined by those 
hard band spectra, we fixed the column density IVh to 
5 x 10 21 cm -2 , which is close to the values obtained by 
Eho et al. (2002) for the SW regions. Note that following 
results have no strong dependence on the column density. 
The standard RMF files version 2006-02-13 were used, 
whereas ARF files were produced by the xissimarfgen 
software version 2006-05-28 with CALDB version 2006- 
05-24 assuming uniform emission from the spectral region. 
As a result, we obtained an acceptable fit fo 2 /degree of 
freedom (dof) = 270.7/291) with the best-fit parameters 
given in Table 2. The best-fit model is shown in Figure 4. 
Note that there is a systematic error of about 15 eV in 
the center energy of the lines. The 7.1 keV line is de- 
tected clearly for the first time. 

Although the model is reproducing the overall spectrum 
quite well, we still see some data excess against the best-fit 
model at ~5.4 keV in Fig. 4. When we added a Gaussian 
with a fixed a (40 eV) to the model, the fit was improved 
with x 2 /dof = 258.5/289. The best-fit center energy and 
the flux of the Gaussian were 5.38 (5.32-5.43) keV and 
2.7 (1.5-3.9)xl0~ 6 photons cm -2 s _1 , respectively. If we 
apply the F-test, the significance level of this additional 
Gaussian is ~ 99.9%. The center energy is consistent with 
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Fig. 2. Count-rate ratio map between the 6.4 keV line and 
the 5.0-6.2 keV band overlaid with contours of the 3.0-6.0 keV 
intensity map. By multiplying 3, a ratio can be converted to 
an equivalent width (keV) of the 6.4 keV line at that position. 



Fig. 3. True color XIS image of the RCW 86 south-west 
region: red represents 0.5-1 .0 keV photons; blue represents 
3.0-6.0 keV photons; and green represents 6.4 keV emission. 


the Ka-line energy of neutral chromium (5.41 keV). 

Since existence of neutral iron at the SW region is sug- 
gested by the Fe-K emission lines, and it may cause ab- 
sorption, we further tried adding Fe K-edge to the model 
with the energy fixed to the center energy of the 6.4 keV 
fine plus 0.712 keV. As a result, the fit was slightly im- 
proved (x 2 /dof = 255.6/288) with the iron column density 
lVp e = 1.7 (0.88-3.5) xlO 18 cm -2 . If we apply the F-test, 
the significance level of existence of Fe K-edge is ~ 93%. 

In order to distinguish between instrumental and intrin- 
sic width, we compared them with those derived from an 
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Fig. 4. Background-subtracted spectra in the hard 
energy band (>3 keV). Data points of XIS 0+2+3 
and XIS 1 are shown with black and red crosses, 
respectively. Solid-lines show the best-fit model 
of a power-law function plus two Gaussians. 


Table 2. Best-fit parameters of the hard X-ray emis- 
sion with a power-law function plus two Gaussians. 


Parameters 


Values 


Power-law 

r 

norm a 

6.4 keV-line 
Center (keV) 
cr b (eV) .... 

Flux c 

7.1 keV-line 
Center (keV) 

<r b (eV) 

Flux 0 

Nr (cm -2 ) ... 
X 2 /dof 


3.17 (3.14-3.19) 
3.5 (3. 3-3. 6) 

6.404 (6.400-6.407) 
61 (56-66) 

5.1 (4.9-5. 3) 

7.08 (7.04-7.11) 
69 (35-101) 

0.53 (0.40-0.69) 

5 x 10 21 (fixed) 
270.7/291 


Notes. Error regions correspond to 90% confidence 


levels. 

Normalization at 1 keV (xlO -2 pho- 
tons keV -1 cm -2 s _1 ). 

This value includes detector energy resolution (see 
text). 

'Total flux (xl0~ 6 photons cm -2 s -1 ) in the line. 


identical analysis done on the Sgr C observation (Table 1) 
which was performed 8 days after the RCW 86 observa- 
tion. The lcr width of 6.4 keV in the Sgr C observation was 
39 (34-44) eV. Since the 6.4 keV line from Sgr C, a molec- 
ular cloud, is expected to have negligible intrinsic width 
compared with the energy resolution of XIS, the width is 
thought to be instrumental energy resolution Therefore, 
the difference of the 6.4 keV line widths between these 
two observations should show the intrinsic width of the 
6.4 keV line from RCW 86, which can be calculated to be 
47 (35-57) eV (99% error region). 




5. Discussion 

The location of the 6.4 keV line has been revealed for 
the first time by the Suzaku observation, and shows that 
the line and the continuum emissions have different origin, 
even though both of them are located at the inner part 
of the shell and are likely to come from ejecta. The solar 
abundance of chromium relative to iron is ~ 0.01 (Anders 
& Grevesse 1989) and if both of Fe and Cr K-a lines are 
excited by electrons (see following discussions), the rela- 
tive strength of these lines should be ~ 0.01. However, the 
observed strength of the Cr Ka line is 0.05 (>0.03; 99% 
region) times that of Fe Ka line. It suggests that the abun- 
dance of chromium was enriched by a supernova explosion, 
especially that of Type la (Iwamoto 1999). Therefore, the 
marginal detection of a Cr-K line also supports the ejecta 
origin of the Fe-K line. 

The Fe-K line features may be fluorescence caused by 
irradiation. Irradiation in this case is the continuum hard 
X-ray emission. We estimate required amount of iron 
by assuming a simple geometry in which the irradiation 
source is a point-like and iron surrounding it is spherically 
symmetrical. Then, the required column density of iron 
Afpe to explain the equivalent width (~ 500 eV) of the 

6.4 keV line is ~ 8 x 10 19 cm -2 . If we assume 4.4 pc (5' 
at 3 kpc) for the radius of the sphere, the mass of iron is 
5.8 x 10 36 g or ~3000 Mq. This is too large for ejecta. 
If we assume solar abundance, the column density of iron 
corresponds to a column density of interstellar medium 
IVh ~ 2 x 10 24 cm -2 . However, no such large molecular 
cloud is reported in the SW region. The column density of 
iron obtained from the Fe K-edge is also discrepant from 
the required value. 

We next discuss the low-ionized plasma scenario for the 

6.4 keV line. From the absence of L-shell lines Borkowski 
et al. (2001a) have pointed out that M-shell electrons have 
to be still bounded. The Fe Kp line, which was detected 
for the first time with Suzaku/XIS, gives more clear con- 
straint on the ionization state of iron. The intensity ratio 
between Fe K a and Kp lines, which is calculated to be 0.10 
(0.074-0.14) (99% error region), is consistent with neu- 
tral iron, and shows that even if iron is ionized the mean 
charge is at most +12 (e.g., Palmeri et al. 2003; Mendoza 
et al. 2004). The center energies of the Fe K a and Kp lines 
are also consistent with this ionization state (e.g., Palmeri 
et al. 2003; Mendoza'et al. 2004). According to Porquet 
et al. (2001), see their Fig. 19, in plasmas of almost all 
temperatures found in SNRs (5 x 10 6 -8 x 10 7 K), iron is 
ionized to a mean charge of +9 — 1-10 when the ionization 
parameter n e t is 10® cm -3 s. If we interpolate the curves 
for n e t= 10 9 and 10 10 cm -3 s, the mean charge of +12 
corresponds to ~ 3 x 10® cm -3 s. Therefore, n e t has to be 
~ 3 x 10 9 cm -3 s or smaller. If we assume the plasma age 
is older than 1000 yr, the electron number density n e in 
the plasma has to be smaller than 0.1 cm -3 . 

On the other hand, there is a required size for the emis- 
sion measure EMp e = fn e nF e dV to explain the strength 
of Fe K a line. If we assume a pure iron non-equilibrium 
plasma (an NEI model; Borkowski et al. 2001b) of kT e = 


1 keV and n e t = 1 x 10 9 cm 3 s, EMp?. has to be 
1.9 x 10 6O d 2 3kpc cm -3 . The elliptical extraction region 
of 12f5 x 9.'2 corresponds to 11x8 pc 2 at 3 kpc dis- 
tance. Assuming a uniform density of a prorate sphere 
with 11 x 11 x 8 pc 3 diameters, we arrive at n e n p e = 
0.006 d^ kpc cm -6 . If we apply the upper limit of the elec- 
tron density, the density of iron np e has to be larger than 
0.06 cm -3 d^ kpc . This density corresponds to an iron mass 


of ~ 8 x 10 34 d% k g or ~40 M 0 , which is too large for 
an ejecta. Therefore, the electron temperature has to be 
higher than 1 keV. If kT e = 2 keV is assumed, both the 
required EMp e and iron mass are about 60 times smaller 
than those in the 1 keV case, and can be acceptable. The 
equivalent width of the Fe Ka line against the contin- 
uum is 530 (480-580) eV. If we assume an NEI plasma 
of kT e = 2 keV and n e t — lx 10 9 cm -3 s, this equiva- 
lent width corresponds to an iron abundance of 1.5. Since 
most of the detected continuum emission is thought to be 
from a different region from the Fe Ka line, this abun- 
dance sets a lower limit. Although the Suzaku observa- 
tion found new aspects on the Fe-K line features and gave 
more concrete limits on the physical parameters, the non- 
equilibrium plasma of Fe-rich ejecta proposed by Bamba 
et al. (2000) and Borkowski et al. (2001a) is still valid for 
the Fe-K line origin. 

Even though the 6.4 keV line can be explained with a 
low-ionized plasma of ejecta, this scenario requires rather 
large mass of iron. Since electrons of ~100 keV energy 
have the smaller cross section for the outer orbital elec- 
trons and the larger cross section for the inner shell elec- 
trons than ~ 10 keV electrons, if we assume harder pop- 
ulation for electrons, the required n e t and n e n F e become 
larger and smaller, respectively, and therefore the required 
amount of iron can be moderated. A non-thermal (supra- 
thermal) spectrum is the hardest case for the electron pop- 
ulation. For example, if we assume electrons in a rather 
hard power-law of index 0.3 with cut-off at 90 keV which 
are assumed for the galactic ridge emission by Valinia et 
al. (2000), the required size of n e npe to explain the 6.4 keV 
line strength is 6 x 10~ 6 d^ kpc cm -6 . This is 1000 times 
smaller than the case of kT e = 1 keV, and therefore the 
lower limit of iron mass is less than 0.1 M 0 . Since elec- 
trons in such hard population would accompany X-ray 
continuum via bremsstrahlung which is harder than the 
detected X-ray continuum (a power- law of photon index 
~3.2), X-ray analysis of the existing HXD data may de- 
tect such hard X-rays, whose spatial correlation to the 

6.4 keV line may be shown by future instruments. 

If only collisional interaction is taking place between 
electrons and ions, kinetic temperature of them are not 
equilibrated with each other at the ionization parameter 
of ~ 10 9 cm -3 s (Itoh 1984), and the ion kinetic tem- 
perature is expected to be 10 times or more larger than 
the electron temperature, or > 10 keV. This high tempera- 
ture suggests a high shock velocity of >3000 km s _1 . This 
value is discrepant with the blast-wave velocity of ~400- 
900 km s~ 1 determined from observations of Balmer- 
dominated shocks (Ghavamian et al. 2001). Plausible 



shocks in the inner part of the shell may have higher 
velocity than the forward shock. The intrinsic width of 
~ 47 eV found for the 6.4 keV line may be evidence of 
this high velocity. The forward shock might have decel- 
erated by colliding with a cavity wall, whereas the inner 
shock might be remaining inside the cavity, or a reflected 
shock after the collision may be moving into the inner re- 
gion. The cavity wall interaction has also been suggested 
by Vink et al. (1997) to explain the overall morphology of 
RCW 86. 

We revealed that the hard X-ray continuum in the inner 
part exists separately from the 6.4 keV line, and there- 
fore confirmed it to be synchrotron emission. In order 
to meet the condition for synchrotron X-ray to be emit- 
ted under the standard diffusive shock acceleration (e.g., 
Drury 1983), the shock velocity has to be larger than 
~2000 km s _1 independent on the magnetic field strength 
(Aharonian & Atoyan 1999). A high velocity shock in the 
inner part of the shell is also required to explain the syn- 
chrotron X-rays. 

An observation with higher energy resolution on the SW 
region is highly encouraged. More strict constraint on the 
iron ionization state can be obtained by measuring the 
difference between center energies of the 6.4 and 7.1 keV 
lines at several eV resolution (e.g., Palmeri et al. 2003; 
Mendoza et al. 2004). Intrinsic width of the 6.AkeV line 
and existence of Cr-K line can also be verified clearly. 

6. Summary 

With the large effective area and low background of 
XIS in the hard X-ray band, we obtained following obser- 
vational results on the hard X-ray continuum and Fe-line 
features from the SW region of RCW 86. 

1. Spatial distribution of the 6.4 keV line is revealed 
for the first time. 

2. Fe K/3 line is detected for the first time. 

3. Strength of Fe K-edge is determined. 

4. We found a sign that Fe Ka line is intrinsically 
broadened for about 50 eV. 

5. Cr-K line is marginally detected. 

All these results are consistent with that the 6.4 keV 
line is emitted from a low- ionized plasma of Fe- rich ejecta. 
We revealed that the hard X-ray continuum is located 
separately from the 6.4 keV line, and therefore confirmed 
it to be synchrotron emission. 

M.U., R.S., H.N., and H.Y. are supported by JSPS 
Research Fellowship for Young Scientists. JPH acknowl- 
edges support by NASA grant No. NNG05GP87G. 
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